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Abstract 

The LHC studies on the MSSM slepton sector have mostly been focused on direct slepton 
Drell-Yan pair production. In this paper, we analyze the case when the sleptons are lighter than 
heavy neutralinos and can appear in the on-shell decay of neutralino states. In particular, we 
have studied the xtx2 associated production, with the consequent decays of xf ~^ ^i^Xi arid 
X2 ~^ ^^Xi via on-shell sleptons. The invariant mass of the lepton pairs, from the neutralino 
decay has a distinctive triangle shape with a sharp kinematic cutoff. We discuss the utilization 
of this triangle shape in distribution to identify the slepton signal. We studied the trilepton 
plus missing Et signal and obtained the effective cross section, a x BR x acceptance, that is 
needed for a 5a discovery as a function of the cutoff mass for the LHC with center of mass energy 
14 TeV and 100 fb~^ integrated luminosity. Our results are model independent such that they 
could be applied to other models with similar decay topology. When applied to the MSSM under 
simple assumptions, it is found that with 100 fb~^ integrated luminosity, a discovery reach in 
the left-handed slepton mass of about 600 GeV could be reached, which extends far beyond the 
slepton mass reach in the usual Drell-Yan studies. 
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I. INTRODUCTION 



While the Large Hadron Colhder (LHC) has great potential in searching for strongly 
interacting particles, its reach in the electroweak sector of new physics scenarios is limited 
due to the suppressed electroweak production cross sections. Although those electroweak 
particles could appear in the cascade decay of heavier colored objects, the discovery reach 
depends strongly on the mass scale of the colored ones. Current LHC searches in the 
Minimal Supersymmetric Standard Model (MSSM) are mostly focused on the direct pair 
production of squarks and gluinos. Null search results from both ATLAS and CMS [l|-l5j 
already exclude the mass of those colored particles up to about 800 GeV. It is likely that 
those colored particles are so heavy that they are out of the reach of the LHC. In this paper, 
we consider the direct production of the electroweak sector of the MSSM. In particular, we 
focus on the LHC reach for the discovery of the sleptons. A complementary study on the 
LHC reach in the neutralino and chargino sector with decoupled sleptons can be found in 
Ref. 0. 

If low energy supersymmetry is realized in the nature, sleptons are likely to be light. 
This happens in the Gauge Mediated Supersymmetry (SUSY) breaking scenarios ^], as 
well as the Anomaly Mediated SUSY breaking scenarios |8| , in which the slepton masses are 
proportional to the electroweak gauge couplings. Even in the minimal Gravity Mediated 
SUSY breaking scenarios (mSUGRA) \^ where all the scalars have a common mass tuq 
at some high energy scale, renormalization group running to low energies typically pushes 
up the squark mass (due to the contributions of strongly interacting gluinos) while the 
sleptons remain light. While the masses of the superpartners of the colored objects have 
already been constrained by current search limits, it is timely to fully explore the discovery 
potential of the LHC for the superpartners of leptons. 

In the i?-parity conserving MSSM with the lightest neutralino Xi being the lightest 
supersymmetric particle (LSP), Y^is a good candidate for Weakly Interacting Massive 
Particle (WIMP) dark matter [lO|. When sleptons are light, the t-channel diagram 
mediated by the exchange of sleptons is important in determining the annihilation cross 
section of Xi dark matter Therefore, discovery of the sleptons is not only a verification 
of low energy supersymmetry in nature; precise measurement of their masses also plays an 
important role in determining the relic density of the neutralino LSP. 

Earlier studies on the slepton discovery potential at the LHC mostly focused on the 
Drell-Yan pair production of slepton pairs, with sleptons directly decaying down to lepton 
and Xi 12l-ll4|. Most of those studies are done either in the mSUGRA framework or for a 
certain set of benchmark points only. The Drell-Yan production cross sections for slepton 
pairs are typically small, suppressed both by the electroweak coupling strength, as well as 
the scalar nature of the sleptons. The LHC reach is very limited: m^^ > 300 GeV and 
mj„ > 200 GeV for the LHC with center of mass 14 TeV and 30 fb~^ integrated luminosity 
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In our study, we focused on an alternative production mechanism for the sleptons via the 
on-shell decay of heavier neutralino and chargino states [l5|. In particular, we considered 
the scenario where Mi < mi < M2 <C fi and studied the pair production of Wino-like 
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X^X2 with the subsequent decay of X2 ~^ ^^X? and Xi ~^ tv.,viL ^J^Xi- The 

colhder signature is trilepton plus missing E^. 

Compared to the traditional searches of slepton Drell-Yan pair production with dilepton 
plus missing Et signature, this channel is advantageous for the following reasons. Firstly, 
the production cross section for Xi'X^' although also at the electroweak strength, is larger 
than slepton pair production due to the fermionic nature of the neutralinos and charginos. 
Secondly, for X2 and Xa being dominantly Wino, the decay branching fraction into left- 
handed sleptons is almost 100%. Thirdly, the SM backgrounds for this trilepton signature 
are also much smaller, with dominant contributions from the leptonic decay of WZj^* ^ as 
well as tt with an extra lepton from heavy flavor decay. 

In addition, the invariant mass of the dileptons from the X2 cascade decay chain, 
m«, has a distinctive triangle sha pe, with a sharp kinematic cutoff, mcut completely 
determined by m^^^ m^o and m^^ [l6|. This triangle feature has been mostly used as 



a precise determination of the slepton mass |16l4l8l|. It is also used in the recent CMS 
opposite sign dilepton searches to enhance the signal acceptance jl]. Although this spectral 
shape can be obvious to the eye, it could easily get washed out when SM backgrounds 
are considered. In our study, we explore the LHC discovery potential for sleptons by 
performing a fit to this triangle spectral shape. We also fit the dilepton invariant mass 
distribution from the dominant SM backgrounds (either containing a Z/'-f* or from tt). 
Since we keep the overall normalization of the backgrounds that contain fitting 
parameter, our treatment allows us to include other backgrounds that contain a Z as well, 
for example, those from other SUSY processes. We obtain the effective cross section, ax 
BR X acceptance, necessary for a 5cr discovery as a function of the cutoff mass. This result 
is model independent and can be applied to any new physics model that gives rise to the 
same decay topology and signature. When applied to the MSSM slepton produced in the 
cascade decay of Wino-like xf X2' obtain the 5a reach in the parameter space of M2 vs. 

A recent analysis by ATLAS on the same sign dilepton plus missing Et signature j2| 
studied x^X2 associated production with decays of Xi and X2 ^ia on-shell slepton in a 
simplified weak gaugino production model. With 1 fb~^ at the 7 TeV LHC, masses of xf 
(X2) up to 200 GeV for m^o = GeV were excluded at 95% C.L., based on event counting. 
Fitting the dilepton invariant mass distribution was also performed in the opposite sign 
dilepton plus missing Et search at CMS [5|. With 0.98 fb~^ integrated luminosity at the 
7 TeV LHC, 95% C.L. upper limits on the cross section times acceptance of about 4 — 30 
fb is obtained for the cutoff mass scale between 20 to 300 GeV. 

The outline of the paper is as follows. In Sec. [TTl we discuss the slepton production 
and decay, focusing on the sleptons produced via neutralino/chargino cascade decay and 
identify their collider signature. In Sec. Illlt we present the triangle spectral shape of the 
dilepton invariant mass distribution and discuss the parameter dependence of the cutoff 
mass. In Sec. HVl we review the current collider bounds on sleptons as well as bounds 
on neutralinos and charginos in the MSSM. In Sec. |Vl we discuss in detail our treatment 
of the signal, as well as the SM backgrounds by fitting to the spectral shape of the m« 
distribution. In Sec. IVTl we present the mo del- independent least effective cross section. 



3 



a X BR X acceptance, that is needed for a 5cr discovery as a function of cutoff mass at the 
LHC with center of mass energy 14 TeV and 100 fb~^ luminosity. In Sec. IVIH we apply our 
study in the MSSM electroweak sector and show the 5cr reach in M2 vs. m^^ parameter 
space. In Sec. IVIIIt we conclude. 



II. MSSM WITH LIGHT GAUGINOS AND SLEPTONS 

We consider the low lying spectrum of the MSSM electroweak sector, which includes 
only neutralinos, charginos, and light sleptons. In our discussion below, we assume the 
canonical case where Mi < M2,n {Mi, M2 and n being the mass parameters for Bino, 
Winos and Higgsinos, respectively) and the lightest neutralino LSP is dominantly Bino- 
type. Scenarios with other mass orderings can be studied similarly. The slepton mass 
spectrum is determined by (M|)2,l, {'Wl'j)^^ and (M|)2,_r,, where each M| is a 3 x 3 matrix, 
representing three generations. In our study, we take the simple assumption that the flavor 
mixing between generations is negligible. The phenomen ology a nd implication of sizable 



flavor mixing in the slepton sector can be found in Ref. [19|, |20|. Furthermore, the left- 
right mixing in the slepton sector is typically proportional to the lepton Yukawa, which 
is small for the first two generations. Therefore, we assume there is no left-right mixing 
for selectrons and smuons and label the mass eigenstates as El and in, for i = e, fi, 
with masses m^^ and "^1^, respectively. For the stau, left-right mixing could be sizable, 
especially with large tan/3. The mass eigenstates are labeled as fi and f"2. There are three 
parameters involved for the stau sector: rrif-^, rrif^ and the left-right mixing angle 6'^. Our 
discussion below applies to the simplified case of selectrons and smuons, although it could 
be adapted to the stau case as well. For sneutrinos, their masses rrij}^ are also determined 
by (M|)|^. Therefore, their masses are related to m^^ with a small splitting introduced by 
electroweak effects: m?^ = + cos 2/3; for the allowed range of tan /3 > 1, < . 

The direct production channels of sleptons are Drell-Yan pair productions of ^l^l, (-l^i, 
vivi and IrIr. The production cross sections are typically small due to both the electroweak 
coupling strength and the scalar nature of the particles. At the LHC with ^/s = 14 TeV, 
the cross sections vary from 0.2 pb to_0.5 fb for ii^i and I'ii'e, from 0.8 pb to 1.5 fb for 
iLi^£, and from 0.08 pb to 0.2 fb for IrEr, for masses of sleptons in the range of 100 to 500 
GeV [13. 

The decay of right-handed sleptons is quite straightforward, proceeding dominantly into 
iXi- In cases when on-shell decay of Ir into higher neutralino states is open and when 
there is a significant Bino component in higher neutralino states (typical for fi ~ Mi), 
^X2 3 4 could also contribute, although the decay branching fraction is almost always 
suppressed. 

The decay of left-handed sleptons depends on the neutralino and chargino spectrum, 
in particular, that of Wino-type neutralino and charginos. Since the decay of sleptons 
into Higgsino-type neutralinos and charginos are typically suppressed, we assume M2 <C /x, 
and thus X2) Xi" are mostly Wino-like. For rn^^, rrii^^ < M2, the branching fractions of 

— ^ ^Xi) ~^ ^iXi are almost 100%. Once m^^, m^;^ > M2, the decays of II iX2^ ^iXii 
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z/£ — )■ z/£X2? ^xf become dominant. The branching fraction is about 10% into Xi? 30% into 
X2, and 60% into Xi- With the subsequent decay of X2 ~^ ^X? and Xi W^*^Xi^ 

left-handed slepton and sneutrino decay would have multi-lepton, multi-jets, plus missing 
Et final states. 

For the Drell-Yan pair production of sleptons with dominant^ direct decay of sleptons 
into Xi, the collider signatures are dilepton plus missing Et for iiii and iniR, single lepton 
plus missing Et for iii^e, and missing Et only for uli^l- The single lepton channel suffers 
from large SM backgrounds, mainly W. The missing Et only signature from Dli^l needs an 
extra jet or lepton from initial or final state radiation, which leads to further suppression of 
signal cross sections. Current collider analyses of slepton Drell-Y an p roduction focus on the 
final states of two isolated energetic leptons plus missing Et 13|, [ij] . The SM backgrounds 



are typically large, dominantly from WW or tt. The LHC reach is very limited: m^^ > 300 
GeV and > 200 GeV for the LHC with center of mass 14 TeV and 30 fb~^ integrated 
luminosity |l4| |. 

In this paper, we explore alternative production channels for sleptons, in particular, ii 
and Ui, via the decay of heavier neutralinos and charginos. The coupling of Higgsinos 
to sleptons are highly suppressed by the small lepton Yukawa coupling, therefore the 
production of sleptons from Higgsino decay is negligible. We assume fi is heavy and 
decouple Higgsinos. Winos, on the other hand, could dominantly decay into ii and z/£ 
once it is kinematically available, since the competing processes of X2 ~^ ^Xi? ^Xi and 
xf — )■ Wxi suffer from small neutralino mixing. The pair production cross sections of 
Wino-like neutralinos/charginos are larger compared to those of sleptons of similar mass. 
In Fig. [H we show the production cross section for Wino-like xfx2 the LHC with 
a/s = 14 TeV. The cross section is about 10 pb for M2 around 100 GeV, which drops 
to about 10 fb for M2 around 700 GeV. In principle, £l,r could also appear in Bino-like 
neutralino decay, in cases of Mi > mi,M2,fi. However, the pair productions of Bino- 
type neutralino with other neutralinos/charginos are typically suppressed due to the small 
neutralino mixing effects. 

For Ml < m^^ , m,>^ < M2 , the lightest chargino Xi" dominantly decays into iui and 

ugiL- With the consequent decay of ui and directly into x?; the branching fraction 
of xf ~^ ^^iXi is almost 100%. X21 the other hand, decays into v^Ui and Ul with 
about the same branching fraction. The former decay leads to vvx^ final states, while the 
latter process has two isolated charged leptons (!,lx\- Considering trilepton plus missing 
Et signatures, the overall branching fraction of xf X2 i^to this final state is about 50%. 
Combining the production cross section of xf X2' left-handed sleptons could be produced in 
the decay products of Wino-like heavier neutralino and charginos states with relatively large 
cross sections compared to the direct Drell-Yan process. Relatively small SM background 
for the trilepton final state and the distinctive triangle spectral shape for m« render this 
channel useful in probing the left-handed sleptons at the LHC. Of course it should always 
be kept in mind that such slepton production is only possible when the slepton masses are 
less than M2. 

It should also be noted that such slepton production via cascade decay of heavier 
neutralinos and charginos only works effectively for the left-handed charged sleptons. For 
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FIG. 1: The left plot shows the cross section for Wino-like X1X2 associated production at the 
LHC with center of mass energy 14 TeV. Here we have decoupled both Higgsinos as well as 
squarks. The right plot shows the Feynman diagram that gives rise to the trilepton plus missing 
Et final states. 



the right handed sleptons, even if it is lighter than M2, the branching fraction is small in 
general compared to the dominant decays of X2 ~^ ^Xi^ X2 ~^ ^Xi since it is suppressed by 
the small Wino-Bino mixing. It only becomes important in the limited parameter regions 
with small M2 — Mi such that the on-shell decay of X2 i'^to the dominant channels are 
forbidden. Therefore, in our study below, we focus on the light left-handed sleptons. 

In cases when /i is lighter: Mi < fi < M2, additional decay modes of heavier Wino 
states into lighter Higgsino states plus Higgses open, which could lead to the suppression 
of the branching fractions of Winos decaying into sleptons. The results that we obtained in 
our study, however, can also be applied to such cases, taking into account the suppressed 
branching fractions. 



III. mu DISTRIBUTION AND TRIANGLE SHAPE 

One distinctive feature of the dileptons from the on-shell cascade decay of X2 (see the 
right plot of Fig. [1]) is that the invariant dilepton mass m« distribution has a triangle shape 
[16], with a cutoff mass determined by the masses of the X12 ^"^^ ^• 




Fig. |2] shows the dependence of mcut on for a given set of (Mi, M2). mcut varies in the 
range of to M2 — Mi. For M2 — Mi > mz, fricnt is larger than mz (indicated by the 
straight red dotted line in Fig. for a large range of m^. This is advantageous since we 
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FIG. 2: The dependence of the dilepton invariant mass m^^ distribution endpoint rricut on the 
slepton mass m^, for (Mi, M2)=(50,400) GeV (sohd curve) and (50,225) GeV (dashed curve). 
Benchmark point I: (Mi, M2, m^-^) = (50, 225, 200) GeV and benchmark point II: (Mi, M2, m^-^) = 
(50, 400, 220) GeV are displayed as dots. Also plotted is the Z boson mass in red dotted line. 



can effectively suppress the dominant SM background from WZ/'y* by imposing a m« cut 
to be above mz- 

This feature in the m« spectral shape is often used as a precise determination of the 



slepton mass [16|. Even in the case of off-shell sleptons in neutralino decay, there have been 



studies in the literature exploring the sensitivity of the m« spectral shape on the slepton 



mass [17|, Il8[. In our study below, we explore how to use this distinctive triangle shape of 
the m££ distribution to identify the slepton signal from the SM backgrounds. This triangle 
shape in m« is not unique to the specified X2 decay in the MSSM; it could appear in many 
new physics model with a similar cascade decay topology that gives rise to two leptons. 
Our analysis is therefore model independent and can be applied to a more general set of 
new physics models. 



IV. CURRENT COLLIDER SEARCH LIMITS 

The current best limits on the slepton masses come from LEP searches for dilepton plus 



missing energy signatures j2l|] with ^/s up to 208 GeV. For a mass splitting between slepton 
and neutralino LSP above 15 GeV and considering only the contribution from right-handed 
sleptons, the mass limits are: mg > 99.6 GeV, > 94.9 GeV and > 85.9 GeV. This 
is conservative, since the production cross section for the left-handed sleptons is higher. 
For stau, it is possible have a large left-right mixing, which could decrease the production 
cross section for the lightest stau pair. A lower limit of rrif > 85.0 GeV can be obtained 
when the production cross section for the lightest stau is minimized. It should be noted 
that the slepton mass limits are obtained with /x = —200 GeV and tan (3 = 1.5, a point at 
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which the neutrahno mass hmit based on the LEP neutrahno and chargino searches is the 
weakest, and the selection cross section is relatively small. The gaugino mass unification 
relation Mi = (5/3) tan^ ^VF^2 is assumed, which is relevant in fixing the masses and field 
content of the neutralinos. Slepton mass limits would change for a non-unified mass relation 
between Mi and M2, since neutralinos appear in both the slepton decay final states, as 
well as participating in the t-channel diagram for selectron production. For selectrons, 
clCr production is also possible via t-channel neutrahno exchange. In the case where the 
C-R ~ Xi mass splitting is small and the usual acoplanar dilepton search is insensitive, a 
single lepton plus missing energy search yields a lower limit on of 73 GeV, independent 
of m^o |22|. For sneutrinos, a mass limit of 45 GeV can be deduced from the invisible Z 



decay width [23[. An indirect mass limit on sneutrinos could be derived from the direct 
search limits on the charged slepton masses. 

Since we consider the production of sleptons from heavier neutrahno decay, we also 
summarize the current status of the neutrahno and chargino sector. Charginos xf can be 
pair produced at LEP via s-channel exchange of Z/'^* or t-channel exchange of Ve-, with 
destructive interference. It decays to ff'Xi via a virtual W or sfermion, or dominantly to 
//' when two body decay is kinematically accessible. In the case of heavy sfermions and 
a mass splitting m ± — m o of at least a few GeV, a robust chargino mass lower limit of 

Al A^l J 

103.5 GeV can be obtained for sneutrino masses larger than 300 GeV 2J], assuming the 
gaugino mass unification relation. For the case of small mass splitting between the lightest 
chargino and neutrahno LSP, limits have been obtained for the degenerate gaugino region 
(Ml ~ M2): m ± > 91.9 GeV for large sneutrino mass, as well as the "deep Higgsino" 



region (|/i| ^ Mi 2): > 92.4 GeV [25|. For lower sfermion masses, the limit is worse 



due to the reduced pair production cross section, as well as the reduction of selection 
efficiency due to the opening of two body decay channels. In particular, there is a so called 
"corridor" region where m^± — mj>^ is small and the lepton from xf ^ ^^i is so soft that 
it can escape detection. Associated production of X1X2 can be studied in cases where the 
chargino search becomes ineffective. Limits on chargino and neutrahno masses for the light 
sfermion case, therefore, depend on the sfermion spectrum. 

For the lightest neutrahno LSP, there is no general mass limit from LEP if the gaugino 
mass unification relation is not imposed. Production via s-channel exchange of Z/'^* could 
be absent for a Bino-like neutrahno, and t-channel production could be negligible for heavy 
selectrons. An indirect mass limit on the neutrahno LSP can be derived from chargino, 
slepton and Higgs boson searches, when gaugino mass (and sfermion mass) unification 
relation is assumed. A lower mass limit of 47 GeV can be obtained at large tan/? (26| . 
while a tighter limit of 50 GeV can be derived in the mSUGRA scenario j27| . 

Trilepton searches at Tevatron Run II jlsl 29 1 study the associated production of xf X2 
with the subsequent decay of Xi ~^ ^^Xi X2 ~^ ^^^~Xi- ^ BR(xfx2 ~^ 3£) is 
bounded to be less than about 0.13 — 0.06 pb (0.5 — 0.1 pb) from D0 (CDF) searches for 
chargino mass in the range of 100 — 180 GeV. For sufficient light sleptons, the leptonic 
decay branching fractions are large and a mass limit on the lightest chargino can be derived 
based on the null search results. A chargino mass limit of 138 GeV is obtained based on 
D0 searches, when the leptonic branching fraction for three body decay is maximized. 
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while no mass limit can be derived in the large mo case 29j. A recent CDF analysis with 



5.8 fb~^ data gives a mass limit on the chargino to be 168 GeV at 95% C.L., for the 
mSUGRA benchmark point mo = 60 GeV, tan/3 = 3 and Aq = (28| . 

Limits on 0" X BR(> 3£) are also obtained based on the recent trilepton search from CMS 
collaboration using 2.1 fb^^ data collected at the LHC with ^/s = 7 TeV j^]- No jet veto 
is imposed and the dominant contribution to the trilepton signal is from gluino cascade 
decay. No limit on the chargino mass can be derived based on the direct pair production 
of XiXi 

Recent analyses by ATLAS on the same sign dilepton plus missing Et [2] studied xfx2 
associated production with the consequent decay of Xi^ and xt "^^^ on-shell slepton. 
Assuming = |(m^o + 'm^±), m^± = m^o, masses of xf (X2) 200 GeV for m^o = 

GeV are excluded at 95% C.L with 1 fb~^ data collected at the 7 TeV LHC. For m^o about 

Xi 

50 GeV, the limit on m + , m^o is weakened to be about 150 GeV. 

CMS performed an analysis on the opposite sign dilepton plus missing Er final states 
by looking for the kinematic edge in the dilepton invariant mass distribution [5|. With 0.98 
fb^^ integrated luminosity at the 7 TeV LHC, 95% C.L. upper limits on the cross section 
times acceptance of about 4 — 30 fb are obtained for the cutoff mass scale between 20 
to 300 GeV, assuming the signal efficiency of the LMl benchmark point: mo = 60 GeV, 
mi/2 = 250 GeV, tan /3 = 10, Aq = and /x > 0. 



V. METHOD 



The dominant standard model backgrounds for the trilepton plus missing Et signal are 
WZ/'-f* and ti with a fake lepton (dominantly from b decay). Note that trileptons from 
heavy flavor bottom and charm decay (produced in association with Z/'j*) could also be a 
significant background lO]- As explained below, the overall normalization of backgrounds 



with a Z peak is a fitting parameter in our analysis; backgrounds containing heavy fiavor 
produced in association with Z/ 7* could be included as well. Therefore, we don't simulate 
such heavy fiavor backgrounds in our analyses. 

Many SUSY models can generate a trilepton signal, for example, x^X2 with X2 ~^ Xi^^*'^ 
and 

xf XiW^*^. For large mass splitting of m^o — m^o > niz, the dilepton invariant 
mass distribution from an on-shell Z looks like that of Standard Model WZ. Such SUSY 
backgrounds, if they exist, are included in our fitting to m« from the Z pole, since the 
overall normalization of the Z contribution is a fitting parameter. For small mass splittings 
m^o — m^o < m^, m« is peaked near m^o — m^o. For such dedicated analysis to 

distinguish such off-shell Z contributions from the triangle spectral shape is necessary. 

The first difficulty in reconstructing the trilepton event is the combinatorial ambiguity 
arising from the presence of a third lepton in the final state. To resolve this issue we use the 
standard technique of same-sign subtraction. We construct our invariant mass histograms 
by including both of the opposite-sign pairs of leptons and then subtract the histogram 
of same-sign dilepton invariant mass to have a good approximation to the histogram 
of invariant masses of the correct pair of opposite-sign leptons, which is otherwise not 
experimentally accessible. While other techniques exist to resolve this ambiguity lH, they 
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all sacrifice statistics for purity, and as our intent is to find the general shape, the statistics 
will generally be of more value for us than the purity. 
We model the dilepton invariant mass distribution by: 

= (/triangle * 9) + fz + Oo/tf- (2) 



dmu 

The first term models the triangle shape of the mu from the signal process: 

= ^^^^^^.^^ , (3) 

where Nsig and rricut are fit parameters for the number of counts in the triangle and the 
mass of the cutoff. To take into account the detector effects, we smear the triangle by a 
convolution of the triangle function with a gaussian with variance shown below: 

(/triangle * ff) = /— / dt 6 ^ ftmng\e{mu - t) . (4) 

y/ Zira j-oo 

The second term models the contribution to from processes that involves a Z peak 
using the Breit-Wigner function: 

A r 

^^^^^^^ ^ 2^(m,,-mo)^ + (r/2)^' 

where A,mo,r are fit parameters for the amplitude, centroid, and width of the Z pole. 
Note that we fit the centroid and the width of the Z pole instead of using the SM values 
to account for the smearing effects introduced by the detector resolution. 

The last term models the contribution to from ti trilepton events, where fu is the 
mu distribution from tt dilepton events taken from Monte Carlo, scaled by a factor ao, a 
parameter to represent the fake rate of tt to give three leptons. This is intended to emulate 
a data-driven background understanding, where the dilepton mass distribution is measured 
in a control region and then used to understand the background in the signal region. We 
choose the dilepton distribution because we expect that the two "wrong" pairs will largely 
cancel each other out in the subtracted distribution. 

Thus, we have in total seven fitting parameters: 

• Number of events in triangle: iVgig 

• Cutoff in triangle distribution: mcut 

• Amount of gaussian smearing of the triangle: a 

• Amplitude of Z peak: A 

• Apparent width of Z peak: F 
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Apparent centroid of Z peak: mo 
Fake rate for tt events: Oq 



We use Madgraph 5 version vO.6.2 [32j and Madevent v4.4.5 7 |33| to generate our signal 
and background events. These events are passed to Pythia v6.4 (34l| to simulate initial state 
radiation, final state radiation, showering and hadronization. Additionally we use PGS4 



35l | with the ATLAS detector card to simulate detector effects. Note that, in producing 
and fitting the dilepton invariant mass, we require only that there be three leptons in the 
event. In particular, there is no requirement of low hadronic activity, which means that 
strong production of particles which later decay through an on-shell slepton can also be 
measured using this technique. We also do not require any missing energy, which allows 
this technique to be applicable in scenarios which do not include an invisible final state 
particle, such as -R-parity violating theories. In general, some minimal requirement will be 
needed to ensure that events can be triggered on (either through a single lepton trigger or 
dilepton trigger), but all of the points we consider have spectra which are not compressed 
enough to have significant loss due to triggering efficiencies. 

For signal generation, we considered the simplified case where the lightest neutralino is 
purely Bino in nature, and the second neutralino is purely Wino, with degenerate charginos 
which are also purely Wino: = Mi and m^± = m^o = M2. This corresponds to taking 
the Higgsino mass parameter /i to be heavy such that the Higgsinos decouple. We assume 
there is no left-right mixing, as well as no flavor mixing between slepton generations. We 
also completely decouple the heavy colored objects. The relevant mass parameters involved 
are Mi, M2 and mi^. 

We simulate the associated production of xfx2 with the consequent decay of xf and 
X2 via left-handed sleptons, as shown in Fig. [1] As discussed earlier in Sec. [Tll since X2 
has equal probability to decay into ii or we obtain trilepton final states 
50% of the time. For simplicity, we only consider trilepton events with £ being either 
an electron or muon. In cases with lepton universality, mg^ = m^^ = m^^, we could 
study the same flavor, opposite sign distribution. For non-degenerate slepton masses 
between generations, multiple triangle shapes appear and the analysis is more complicated, 
though with no flavor mixing they appear in different channels. Note that in the realistic 
case when the neutralino and chargino mass eigenstates are not the pure gauginos, the 
corresponding branching fractions for the decay into Cl and need to be considered. For 
the backgrounds, we generate the SM WZ /Y trilepton events, as well as tt with both tops 
decaying leptonically. 

We construct our signal and background histograms of the dilepton invariant mass by 
using our Monte Carlo data as probability distributions from which we select events. For 
the signal, we draw from the opposite-sign lepton events 2Nsig times and draw from the 
same-sign lepton events A^sig times and construct the difference between the opposite-sign 
and same-sign distributions. A^^sig is the number of events that we are including in the 
pseudoexperiment we are currently generating. For a given luminosity, A^sig = C x cTgig x 
BRx acceptance. Similarly, we can draw N^zj-y* events from the SM WZ/'-j* background 
and construct the corresponding m^^ distribution. 
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FIG. 3: mu distribution for MSSM trilepton signal and dominant SM backgrounds at the 14 TeV 
LHC with an integrated luminosity of 100 fb~^. The relevant MSSM parameters are chosen to 
be benchmark point I: {Mi,M2,m^J = (50,225,200) GeV with mcut =100 GeV (left plot), and 
benchmark point II: (Mi,M2,m^-^) = (50,400,220) GeV with rricnt =325 GeV (right plot). Also 
shown as black curves are the best fit distribution using Eq. ([2]). 

Purely leptonic ti decay could lead to trilepton events with a third faked lepton from 
b decay. In the pseudoexperiment we are currently generating, we estimate the expected 
trilepton events Nu using the fake rate estimated from PGS simulation: Nu = C x au x 
BR X fake rate x acceptance, where the fake rate is about 4.3 x 10~^. 

In a trilepton ti sample, the opposite sign dilepton mu distribution originating from 
W^W~ decay is the same as the dilepton distribution from dilepton tt events. While 
opposite-sign/same-sign dilepton rriei distributions with one lepton from W and the other 
from near/far b jets have a very different distribution in principle, due to poor Monte- 
Carlo statistics on the three lepton ti events, we draw from the opposite-sign dilepton 
distribution of the dilepton ti event samples N^t times to simulate the expected final 
distribution of m« from W pair decay. We take the simplified assumption that the same 
sign and opposite sign mu distributions from the Wb combination are very similar and 
build two random histograms from the same-sign dilepton distribution of the trilepton ti 
sample with N^^ entries each to simulate the wrong pair of opposite-sign leptons and the 
same-sign pair of leptons, respectively. We then combine these histograms appropriately 
to get our pseudoexperiment result histogram for the ti background. While this is not 
a perfect description of the background in question due to physical differences between 
the distributions of the incorrect pair of opposite-sign leptons and the pair of same-sign 
leptons, it is approximately valid and a more accurate data-driven background method 
than the one we use will be able to include those differences without great difficulty, so the 
sensitivities we find using this technique should be valid. 

In Fig. [31 we show the m« distribution for the MSSM trilepton signal and dominant SM 
backgrounds at the 14 TeV LHC with integrated luminosity of 100 fb~^, for two benchmark 
points: I, (Mi,M2,m^- ) = (50,225,200) GeV (left plot); II, (Mi,M2,m^- ) = (50,400,220) 
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A 


r 

(GeV) 


mo 
(GeV) 


ao 
(xlO-3) 


(GeV) 




a 

(GeV) 




Benchmark 
Point I 


input 
fitted 


3.981 X 10^ 
3.299 X 10^ 


3.121 
(3.121) 


90.158 
(90.158) 


4.280 
5.420 


99.8 
99.5 


1.85 X 10"^ 
1.83 X 10"^ 


0.980 


1.54 


Benchmark 
Point II 


input 
fitted 


3.981 X 10^ 
3.347 X 10^ 


3.121 
2.914 


90.158 
90.041 


4.280 
4.412 


325.3 
321.9 


3.17 X 10^ 
3.14 X 10^ 


2.319 


0.94 



TABLE I: Input and final fitting parameters to Eq. ([2]) for two benchmark points: I, 
{Mi,M2,m^J = (50,225,200) GeV; II, {Mi,M2,m^J = (50,400,220) GeV, for the 14 TeV 
LHC with £ = 100 fb^i. The cutoff for benchmark I is less than 150 GeV so we fix the width 
and centroid of the Z peak with the values obtained from the m^i distribution of the SM WZ/'j* 
backgrounds based on Monte- Carlo simulation. 



GeV (right plot). The corresponding triangle cutoff masses are 100 GeV and 325 GeV, 
respectively. For rricut near mz (as in Benchmark Point I), the triangle distribution of mn 
from the signal process is buried under the SM Z pole, making the identification of such 
triangle features much more difficult. For mcnt far above mz (as in Benchmark Point II), 
the sharp cutoff feature in m« distribution can be easily identified from the SM background, 
as shown in the right plot of Fig. |3l 

For a given value of integrated luminosity, we fit the mu distributions built from the 
signal and background events generated as described above, using formulae given in Eqs. ([2]) 
— dSJ with seven fitting parameters. When the cutoff falls near the Z pole, we find 
that the fitting routine has too much freedom and often fails to identify the cutoff in the 
proper location. This failure is due to the similarity of the sharp feature of the cutoff and 
the edge of the Z pole. This results in a large degeneracy in the signal and Z pole fit 
parameters. Therefore, we fix the Z width, F to be 3.121 GeV and centroid position, mo 
to be 90.158 GeV when the cutoff is less than 150 GeV. These values are obtained from 
the m« distribution of the SM WZ/'y* backgrounds for events passed through the PGS 
detector simulation. 

We perform a fit using the MINUIT fitting routine [1^. We fit from 20 GeV to 
("^cut+100 GeV), with a fixed binning scheme of 3 GeV/bin. We find a good fit to the 
data with x^/dof ~ 1 for all cutoff masses. For high cutoff masses, x^/dof is slightly lower 
due to the fact that our fit spans more bins. Conversely, we find a slightly higher x^/dof 
at lower cutoff mass. The data driven approach to fitting the background works very well 
providing a best fit x^/dof of about 1 when fitting only the background. For high cutoff 
masses, we also consider fits from 125 GeV to (rricut+lOO GeV). This to ensure that our 
fits for the cutoff at high mass are not being influenced too greatly by fitting the Z at low 
mass. As an illustration, we show in Table [T] the input parameters for the simulation, as 
well as the fitting parameters for the two benchmark points. 
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FIG. 4: Effective trilepton cross section, a x BR x acceptance, at the 14 TeV LHC required for 
a 5(T detection of the cutoff feature with an integrated luminosity of 100 fb~^. For the left plot, 
the cross section reach is obtained with a fitting for mu in the range of 20 GeV to (m-cut+lOO 
GeV). For the right plot, a fitting range of 125 GeV to (mcut+100 GeV) is used. 

VI. RESULTS 

Using the fitting strategy described above and assuming only the SM backgrounds for 
trilepton plus missing Et signals, for a given cutoff mass mcut, we fit for the cutoff feature, 
marginalizing over the other fit parameters. This allows us to determine the number of 
events required for a detection of the cutoff feature with Sex-level confidence. We consider 
cutoff masses ranging from 75 to 650 GeV stepping by 25 GeV, with finer stepping around 
the Z pole. As explained before, for cutoff masses below 150 GeV, we perform a five- 
parameter fitting with the Z width and centroid fixed. For cutoff masses above 150 GeV, 
we perform a seven-parameter fitting, allowing the fitting parameters of the Z pole to vary 
freely. 

In the left plot of Fig. H] we plot the required effective trilepton cross section, a x BR x 
acceptance, as a function of cutoff mass for 5a discovery at the 14 TeV LHC with an 
integrated luminosity of 100 fb~^. The maximum near 100 GeV is due to the difficulty of 
detecting the cutoff feature near the Z pole. There is a significant parameter degeneracy 
between the amplitude of the Z and the number of counts in the triangle. The effect also 
creates a significant scatter amongst the data points. At cutoff masses greater than 150 
GeV, the required cross section for detection decreases with increasing cutoff mass. The 
scatter at high cutoff mass is due to statistical fluctuations in the Monte Carlo during the 
fitting process. The quality of the fit, however, depends on the range of m^i that is used 
in the fitting. 
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In the right plot of Fig. HI we show the 5a reach in the effective cross section similar 
to the left plot, but only fitting the m« distribution above 125 GeV rather than above 
20 GeV, such that the background contributions are significantly less relevant to the fit. 
Comparing to the left plot, in which a blind fitting is performed without prior knowledge 
of possible range of mcut, although the resulting lower limit on the effective cross section 
is very similar in both cases, the fit with m« > 125 GeV above the Z pole is more robust 
since it greatly reduces the dependence on the precise knowledge of the backgrounds, in 
particular, those containing a Z pole. 

Note that the above result is obtained for an integrated luminosity of 100 fb^^. There is 
no simple scaling behavior of the required effective cross section for a different luminosity 
because this fitting technique does not have the simple statistical behavior of a counting 
experiment. In order to understand the sensitivity of this technique at significantly different 
luminosities it is necessary to generate a new set of pseudoexperiments and analyze them 
as explained above. 



VII. SLEPTON REACH IN MSSM 

The most straightforward application of this search is to a system which gives the 
maximal branching ratio for the slepton decay of the second neutralino. We therefore 
consider the simplified MSSM scenario that is described earlier, in which Xi is purely Bino 
and Xi^) X2 purely Winos with no left-right mixing in the slepton sector and squarks 
decoupled. Note that in this simplified scenario, such a trilepton search channel is only 
sensitive to the intermediate ii since decays to iji are forbidden due to the absence of 
couplings. 

Wino pair production x^X2 is completely controlled by gauge couplings in this 
approximation, which allows us to make robust predictions for the cross section of this 
process. With our fitting strategy, for a given mass parameter set with known production 
cross section, we allow the luminosity to shift until the triangle spectral shape in m« 
distribution is detectable at the 5a level, as defined above. The previous limits from LEP 
II constrain the left-handed sleptons to be heavier than about 100 GeV, and so we consider 
only sleptons which are 100 GeV or heavier. While charginos are more weakly constrained 
in general than the left-handed slepton, they must be heavier than the left-handed slepton 
in order to be within this framework. We scan m^^ and M2 in the range of 100 — 700 GeV, 
with a step size of 25 GeV. The luminosities necessary to discover the slepton in these 
decays as a function of the masses are shown in Fig. [51 

In the left plot of Fig. |5l we assume that only one lepton flavor (selectron in our analysis) 
is accessible in the cascade decay of X2- For the 14 TeV LHC with 100 fb~^ luminosity, 
the mass reach for the ii can be reached up to 600 GeV at 5a level. In the right plot, 
we assume the lepton universality condition such that mg^ = m^^ = m^^. Considering 
only trilepton events with electrons and muons, namely eee, /x/x/i, e+e~/i^ and /i"'"/i~e^; 
additional branching ratio suppression factors need to be applied. The reach for the left- 
handed slepton mass is reduced in this case, about 500 GeV for fri^^. 



For comparison, Drell-Yan slepton searches were considered in Ref. [lj|, with the authors 
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FIG. 5: Luminosity required at the 14 TeV LHC for a 5a detection of the cutoff feature above 
the SM backgrounds. In the left plot, we assume only the left-handed selectron is accessible in 
the cascade decay of X2- the right plot, we assume that all slepton masses are degenerate. The 
cross section reaches are obtained considering the final states with electrons and muons only. 

concluding that flavor-diagonal sleptons with masses less than 350 — 400 GeV can be 
discovered by CMS using 100 fb"^ of data. Thus, our technique indicates that, in the 
pure gaugino limit, the reach for left-handed sleptons will be significantly enhanced by 
looking for the characteristic cutoff shape of m« in the decays of the gauginos. Note that 
right-handed sleptons are not subject to these results because they are a singlet under 
SU(2)i. 



VIII. CONCLUSIONS AND DISCUSSION 

In this work, we have studied the LHC discovery potential for sleptons, which are 
produced via the on-shell decay of the heavier neutralino and chargino states. In particular, 
we have studied the xfx2 associated production, with the consequent decays of xf ~^ ^dXi 
and X2 ~^ ^^Xi via on-shell sleptons. Comparing to the conventional slepton searches 
through Drell-Yan production and dilepton plus missing Et final states, this trilepton 
plus missing Et channel has the advantage of larger production cross sections and less 
SM backgrounds. In addition, the invariant mass distribution of the dilepton pair from 
X2 decay has a distinctive triangle shape, which can be utilized to select out the signals 
from the dominant SM background of W Z/'-y* and tt fakes. We performed a fit to the m« 
distribution of both the signal and the backgrounds. For the LHC with 14 TeV center 
of mass energy and 100 fb~^ integrated luminosity, we obtained the model-independent 
lower bounds on the effective signal cross section, a x BR x acceptance, as a function of 
triangle cutoff mass, mcut, at 5cr significance level. Applying this result to the MSSM in 
the parameter space of Mi < m^^ < M2 <^ /i, we found that the mass reach for the 
can be up to 600 GeV at 5cr level at the 14 TeV LHC with C =100 fb~^, when there is 
only one slepton generation (selectron in our study) lighter than Winos. For a degenerate 
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slepton spectrum with mg^ = uifX^^ — mf^, and considering final states with electrons and 
muons only, the reach is slightly worse due to the suppression of the branching fractions. 

Comparing to earlier studies on the LHC reach for sleptons from Drell-Yan production, 
the reach for left-handed sleptons via Wino decay is greatly enhanced. On the other 
hand, it should be noted that our study works most effectively for the left-handed slepton, 
since the decay fraction of heavier neutralino/chargino states to right-handed sleptons is 
typically suppressed in most of parameter space, cither by the small Bino-Wino/Higgsino 
mixing, or the lepton Yukawa couplings. The right-handed slepton could appear in Bino- 
like neutrahno decay, if it is not the LSP. The associated production cross sections for Bino 
with other neutralino/chargino states, however, are suppressed in general. Our study does 
not apply to the parameter region of m^^ > M2 since the on-shell decay of Winos into 
left-handed sleptons is forbidden by kinematics. Therefore, for left-handed sleptons with 
mass heavier than M2 or for right-handed sleptons, the usual Drell-Yan production is still 
the dominant production mode. 

It should also be emphasized that the results we obtained for the lower bounds on the 
effective signal cross section, a x BR x acceptance, as a function of triangle cutoff mass 
scale mcut is model independent since it can be applied to other new physics models that 
give rise to the same cascade decay topology and final states of trilepton plus missing Et- 
For new physics models with a given parameter set, we can obtain m^ut as well as the 
production cross sections, branching fractions into the trilepton final states, and signal 
acceptance via detector simulation. Comparing it with the lower bounds we obtained, we 
can derive the LHC reach in the parameter space for such new physics models. 

More work is needed to fully explore the LHC reach for the slepton sector. Our results 
on the left-handed slepton reach are obtained under the simple assumption that Xi is a pure 
Bino state, X2>xf pure Wino states, and the heavier Higgsino states are completely 
decoupled. In addition, we studied only final states including electrons and muons. The 
analysis strategy in our study can be applied to the general MSSM framework with the 
mixing of gaugino and Higgsino states, as well as three lepton flavors (and possible left- 
right mixing in the stau case). The corresponding branching fraction into trilepton final 
states needs to be taken into account in such general cases. Trilepton final states with 
all possible flavor combinations could also be studied, although the fltting to the triangle 
shapes might be more complicated when the slepton masses are not degenerate. 

In our study, we performed the triangle shape fitting to the trilepton plus missing Ej^ 
final states. The triangle shape in ma distribution arises from the X2 cascade decay chain. 
Therefore, it appears in any process that contains such a heavy neutralino cascade decay. 
The fitting strategy could be apphed to final states containing two opposite sign same 
fiavor dileptons. For example, charginos in x^X2 production could decay to jets instead 
of leptons; or we could consider productions originated from gluinos or squarks, with the 
cascade decay of gluino or quark containing a X2- SM backgrounds for dilepton plus jets 
plus missing Et signature, of course, are very different from the trilepton plus missing Et 
signal that we considered in our study. 

We could also consider Wino type xtXi production with xf ~^ ^dx^- The final state 
of dilepton plus missing Et is similar to the conventional slepton study of Drell-Yan pair 
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production oi lili with £l — )■ Although this channel is not as powerful as the trilepton 
plus missing Et study that we explored in this paper, it would still have better reach 
comparing to the Drell-Yan process given the larger production cross sections. 

As we mentioned earlier, our analyses can not be applied to the cases when the sleptons 
are heavier than M2, when the Drell-Yan is the dominant production mode. Previous LHC 
analyses on the Drell-Yan channel focused on the dilepton plus missing Et final states. 
When left-handed sleptons are heavier than Wino-like neutralino/charginos, the branching 
fractions of heavier sleptons into Wino-like states are sizable given its SU(2) coupling 
strength. Considering the consequent decay of the Wino-like neutralino/ charginos, multiple 
leptons (up to six) + jets + missing Et final states could appear, which provides additional 
channels for the left-handed slepton discovery at the LHC. For the right-handed sleptons, 
however, decays into lepton plus Bino-like x? LSP are still dominant. 

Light sleptons could also contribute sizably to low energy processes, for example, parity- 



violating electron scattering, leptonic Pion and Kaon decays, etc. [37[. Given the recent 
progress on both the theoretical and experimental studies in those low energy precision 
measurements, they have reached a sensitivity which is now able to probe new physics 
beyond the SM. LHC studies on the slepton sector will be complementary to the indirect 
probes provided by these precision measurements. 
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